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PRINCIPAL  NOTATION 
Acceleration  vector 
Contravarlant  componenta  of  a 

Orthogonal  surface-imbedded  vectors  such  that  £.  Is  along  the 
constant  f ractlon-of-chord  lines 

General  coordinate  vectors,  vector  components 

Orthogonal  surface-imbedded  vectors  such  that  c^^  Is  along  the 
tangent  to  the  surface  In  the  chordwlse  direction 

Alternating  unit  tensor 

Unit  vectors  In  r  and  6  directions  on  blade 

Unit  vectors  In  helical  coordinate  system 

Combinations  of  geometric  variables 

Metric  tensor 

Associated  metric  tensor 

Cofactor  of  g^j 

Total  rake  (see  Figure  3) 

Chrlstoffel  symbol  of  the  first  kind  (see  Equation  13) 

Unit  normal  to  blade  surface  (see  Figure  4) 

Velocity  vector  and  components 
Position  vector- 
Radial  coordinate 

Position  vector  of  point  on  blade  surface 
Far  field  flow  velocity  (ahead  of  propeller) 

Orthogonal  surface  coordinates  (see  Figure  4) 

Body  force  component 
Coordinate  axes 

Radial  and  chordwlse  parameterp  defining  locations  on  tlie  mean 
blade  surface 
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Conblnatlons  of  geonetrlc  variables 

Angular  velocity  of  rotating  blade 

Vortlclty  vector 

Absolute  viscosity 

Kinematic  viscosity 

Helix  pitch  angle  of  reference  blade 

Skew  angle  (see  Figure  3) 

Chrlstoffel  symbol  of  second  kind  (see  Equation  14) 

Generalized  Kronecker  Delta 

Coordinates  of  blade-oriented  coordinate  system  (see  Reference  1 
and  Figure  3) 
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ABSTRACT 


\  , 

^Specific  expressions  for  a  general  tensor  description  of 
the  full  field  equations  for  flow  about  a  propeller  blade  are 
given.  Two  orthogonal  blade-surface  coordinate  systems  are 
considered:  one  consisting  of  the  normal  and  constant  radius 
lines  and  one  consisting  of  the  normal  and  constant  fraction- 
of-chord  lines.  The  third  direction  is  obtained  by  the  cross- 
product  of  vectors  in  the  two  given  directions.  It  is  believed 
that  the  first  coordinate  system  is  appropriate  to  describe  the 
flow  over  the  majority  of  the  blade,  and  the  second  coordinate 
system  may  be  better  suited  to  describe  the  flow  at  the  leading 
edge.  It  is  recommended  that  the  boundary-layer  limit  of  the 
field  equations  be  determined  and  numerical  evaluations  be  con¬ 
ducted  to  explore  the  merits  of  the  two  coordinate  systems. 

/\ 
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INTRODUCTION 

The  various  performance  properties  of  propellers  have  generated  a  significant 

number  of  geometries  for  different  purposes.  It  has  been  generally  accepted  that 

a  blade-fixed  coordinate  system  is  most  useful  for  general  application.  Prom  this 

conclusion,  one  is,  of  course,  faced  with  the  necessity  of  adding  new  terms  to  the 

1* 

equations  of  motion  dictated  by  the  rotating  coordinate  system. 

The  equations  of  motion  in  their  general  form  have  been  considered  by  many  and 

12  3 

have  been  written  down  in  a  general  way.  *  '  Depending  upon  the  problem  to  be 
solved,  the  simplifying  assumptions  made,  and  the  technique  of  solution,  one  may 
find  some  coordinate  systems  have  advantages  over  others.  It  is  not  the  purpose 
here  to  pass  Judgement  on  any  coordinate  system  but  rather  to  suggest  an  approach 
to  some  hydrodynamic  (fluid  mechanical)  problems  and  offer  the  details  of  appro¬ 
priate  coordinate  systems.  In  at  least  one  case,  a  specific  coordinate  system  is 
offered  that  has  not  appeared  in  the  literature  Involving  propeller  theory  and 
that  is  believed  useful  to  described  boundary-layer  flow  on  the  blades.  There 
are  areas  for  which  the  suggested  coordinates  may  need  significant  modification. 


*  A  complete  listing  of  references  is  given  on  page  37 


such  as  problems  Involving  cavity  flows.  Certainly,  the  boundary  conditions  and 
flow  surface  boundaries  will  need  attention.  In  some  cases,  these  may  Involve  an 
Iteration  to  obtain  a  proper  boundary.  Whatever  Is  needed  can  at  least  be  Incorpo¬ 
rated  Into  the  current  coordinate  systems;  furthermore,  the  problems  themselves  may 
suggest  new  geometries  that  may  be  more  useful. 

The  notation  used  In  the  report  Is  that  of  generalized  tensors.  The  basic 
premise  Is  that  one  can  set  down  the  momentum,  continuity,  and  energy  equations  and 
any  equations  of  state  In  this  general  descriptive  manner.  Once  this  Is  accom¬ 
plished,  any  convenient  geometry  can  be  Introduced.  The  selection  can  be  made  (and 
often  Is)  to  closely  accommodate  physical  conditions,  to  permit  partial  uncoupling 
of  the  equations  (or  at  least  some  simplification),  or  to  describe  boundary  con¬ 
ditions  rather  simply.  The  coordinate  system  so  chosen  may  not  always  be  an  orthog- 
vjnil  one,  but  in  certain  applications  may  be  more  convenient  and  perhaps  more 
economical  to  solve  numerically.  In  the  following,  the  basic  fluid-dynamics 
equations  are  Introduced  in  this  generalized  tensor  form  without  a  great  deal  of 
development.  Details  of  some  of  the  requisites  are  presented  In  Appendix  A,  and 
the  reader  Is  referred  to  Reference  6  for  other  details.  Covariant  and  contravarl- 
ant  forms  of  the  equations  are  Introduced.  Some  emphasis  Is  placed  on  the  contra- 
variant  form  because  most  people  feel  more  comfortable  and  have  a  better  physical 
feel  In  that  form  due  to  the  close  relationship  of  contravarlant  components  with 
physical  components.  Reference  is  made  in  the  present  report  to  advantages  and 
disadvantages  In  particular  cases  of  both  systems. 

The  lengthy  mathematical  manipulations  are  not  Included  In  the  text;  however. 

In  some  cases.  Intermediate  steps  are  Included  to  make  the  procedure  more  easy  to 
follow  for  those  not  familiar  with  the  tensor  notation.  The  tensor  notation  and 
manipulation  are  believed  to  have  some  advantages  over  the  usual  vector  notation  In 
two  ways: 

1.  One  knows  precisely  how  to  expand  a  term  or  operator  as  soon  as  the 
geometry  Is  selected. 

2.  As  often  happens,  higher  rank  tensors,  (e.g.,  stresses)  and  their  deriva¬ 
tives  require  no  more  expertise  to  handle  than  first  rank  tensors  (vectors),  and 
the  concept  of  dyadlcs,  which  Is  unfamiliar  territory  to  many,  need  not  be  employed 

The  present  work  Is  Intended  to  provide  a  specification  of  geometry  Immedi¬ 
ately  useful  for  some  hydrodynamic  propeller  problems.  In  particular,  the  full 
field  equations  are  given  for  flow  about  a  propeller  blade.  These  field  equations 


can  be  evaluated  to  determine  which  terms  should  be  retained  for  boundary-layer  flow 
on  the  blade  surface*  It  Is  hoped  that  the  manipulations  described  will  provide  an 
easy  transition  to  other  areas  of  endeavor  as  well* 

ANALYSIS 

TENSOR  DESCRIPTION  AND  PHYSICAL  EQUATIONS 

While  It  Is  not  necessary  to  restrict  the  equations  of  hydrodynamics  to  Incom¬ 
pressible  fluids,  this  approximation  Is  employed  here  simply  to  address  the  specific 
purpose  Involved  In  this  report,  namely,  problems  associated  with  hydrodynamics  and 
even  more  specifically,  the  hydrodynamics  of  propellers*  The  essential  equations 
are  given  below  In  contravarlant  form: 

Continuity  Equation: 

q*.  =  0 

^  '  I  alVgq'l 

"Viax*  (1) 

Equation  of  Motion: 


a+  ^xR  +2nx^-^nx[nxRl=F-^Vp 


(2) 


where  F  is  the  system  of  forces  such  as  body  forces  or  viscous  forces*  The  terms 
containing  Q  are  due  to  the  rotating  coordinate  system  chosen.  Because  I?  Is  con¬ 
sidered  to  be  constant,  the  time  derivative  term  Is  zero.  The  2l2cq  term  Is  the 

Coriolis  acceleration  due  to  the  accelerating  (rotating)  system  and  the  last  term 

4 

is  the  centripetal  acceleration*  The  acceleration  measured  in  the  chosen  coordi¬ 
nate  system  is  a,  with  the  following  components: 


(3) 


*The  summation  convention  here  Is  sum  on  repeated  sub  and  superscripts. 


3 


where  Is  the  Intrinsic  (substantial)  derivative*  For  steady  state,  this 

reduces  to 


a*  = 


A) 


i 


If  one  looks  at  the  viscous  desc)riptlon  of  Reference  5,  it  is  noted  that  the 
viscous  stresses  are  not  changed  because  of  the  rotating  system  of  coordinates* 
Physically,  this  makes  sense  If  one  recalls  that  all  viscous  stresses  Involve 
velocity  gradients  and  not  the  velocities  themselves* 

Additional  terms  do  arise  if  one  considers  turbulent  motion  and  the  type  of 
term  associated  with  perturbations  about  a  mean  flow*  In  the  case  of  turbulent 
flow,  an  additional  term  will  appear  an  additional  term  due  to  Reynolds*  stresses. 
Secondary  motions  are  not  excluded  In  Equation  (2)  and  may  play  an  Important  part 
In  some  flows,  although  they  normally  do  not  become  Important  except  In  the  last 
stages  of  an  axial  flow  compressor  or  In  the  root  regions  of  a  propeller. 

One  significant  term  that  needs  some  discussion  is  the  vortlclty  vector* 


Vxq  =  5 


(5) 


If  there  were  no  vortlclty  far  ahead  of  the  propeller.  It  would  be  necessary  that  U 
(the  flow  velocity)  be  constant  or  a  function  of  time  only.  If  It  were  not,  if,  for 
example,  U  =  U  (r)i  one  could  not  even  suggest  that  the  far  field  Is  a  potential 
field*  This  would  cause  no  great  problem  In  the  previous  equations  because  there 
is  little  to  be  gained  by  such  an  assumption  If  Indeed  one  wants  to  consider  the 
viscous  terms* 

For  general  flow  of  the  type  described  here,  the  viscous  component  of  F  for  an 
Incompressible  fluid  Is  accordingly 

Ft  =  M9"'q,qK 


The  remaining  terms  to  be  put  In  tensor  form  are  2  flxq  andfixIflxR),  where  R  Is  the 
position  vector  In  the  rotating  system  of  coordinates* 


<*>  \  1% 


Figure  1  -  Centripetal  Acceleration 

Hence,  Is  the  centripetal  acceleration  In  the  negative  r  direction 

vtovard  the  =  x  axis  of  rotation).  The  term  i?xq  Is  found  from  the  following 

1  liTiire : 


Figure  2  -  Coriolis  Acceleration 

The  components  of  this  vector  must  be  determined  In  a  particular  coordinate  system. 
Q  must  be  specified  In  the  particular  system  also.  For  laade-orlented  coordinates 

[IJ: 

ii  =  nsin0pe,+ Hcos^pCj  (7) 

and  for  both  cartesian  and  cylindrical  polar  coordinates: 

n=-nr  (8) 


In  expanded  form  (which  one  must  use  for  solution,  the  equations  of  notion  are 

53'  +  <i'(-|^+  ni's')  +2VS9".,„n'<i‘+«Mnin''*'  = 

-  -ro'i^  d*q'  .  pi  dq'  .  I'i  dq*  .-fdq* 


where  X  are  the  contravarlant  components  of  the  body  forces  (e.g«,  buoyancy). 
The  equations  of  motion  can  be  written  In  covariant  form  as: 


+  =  -  l|£  +wg«'‘q,,,k  +  X,  -2e„fcn'q'«-6|Jn'n,Rfc 


q,  ,^=_i!3j - r°i3j  \ 

axiax**  '•  dx^  "‘dx* 


The  above  equations  have  been  written  for  a  right-handed  coordinate  system*  A  left- 
handed  system  will  change  some  of  the  signs.  The  argument  for  considering  the 
contravarlant  form  of  the  differential  equations  (9a)  as  opposed  to  the  covariant 
form  (9b)  Is  that  only  contravarlant  components  of  the  velocity  are  Involved  In  the 
equations  (9a)  while  one  has  a  mixed  system  in  covariant  and  contravarlant  forms 
(9b).  The  many  symbols  Introduced  above  are  detailed  below; 

e^jj^  Is  the  alternating  unit  tensor  and  has  the  definition 


'ijk  ■  " 


i=jorj=korl»kori=j=k 


•*1  Ifij^k^l  and  IJk  are  cyclic  as  1231 

="-l  l/jfk^l  and  ijk  are  not  as  acyclic  3213 


I 
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X  Is  the  usual  Kronecker  delta: 

■  0  m  #  n 

-  1  Di  »  n 

5^^  has  the  definition 

Opq 


-  1 

1 

1* 

1 

and 

1  -  p 

and 

1  -  q 

-  -1 

1 

1 

and 

1  -  q 

and 

1  -  p 

»  0  for  all  other  cases. 

1 1k 

e  Is  the  same  as  with  relative  weight  +1  Instead  of  -1  (for 

g.  .  Is  the  metric  tensor 
Ij 

Is  the  associated  metric  tensor 


r.  are  the  Christoff el  symbols  of  the  second  kind.  They  are  not  tensors, 

f  1 1 

and  for  this  reason  many  people  prefer  the  notation  Instead  of  the  Princeton 

notation  [6]. 


Physical  vector  quantities  are  related  to  the  forms  as  follows.  The  terms 

v/^  ql(not  summed)  are  the  physical  components  of  q  that  are  the  edges  of  the 
^11^  _  . 

1  iiralleleplped,  of  which  q  Is  the  diagonal.  The  quantity  y Vgn  (not  summed)  is 

the  length  of  an  orthogonal  projection  of  q  on  the  tangent  to  the  coordinates  at  a 

point.  Because  •  ^rs**  ’  obtain  the  covariant  form  by  first  defining  the 

contravarlant  form 


k  — 

where  Q  are  the  physical  components  of  Q,  and  then  deriving  q^^  via  the  above 
equation. 

At  this  point  one  needs  to  develop  the  values  of  the  contravarlant  components 
of  specified  vectors  plus  the  Chrlstoffel  symbols  Indicated.  One  now  must  specify 
the  geometry.  We  shall  carry  out  some  of  the  details  that  we  have  been  indicating 
for  the  general  fluid  flow  problem. 

*See  Section  41  of  Sokolnlkoff  [6]  for  detailed  explanation. 
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COORDINATE  GEOMETRY 

The  (Y  Cartesian)  system  and  the  X^:(x,r,0)  system  are  well  described  In 
standard  texts,  e<g«,  reference  [6].  The  details  are  contained  In  Appendix  B. 
However,  propeller  blades  and  the  helical  nature  of  the  associated  fluid  motion 
require  somewhat  unique  coordinate  systems*  One  such  system  Is  shown  In  Figure  3« 


Figure  3  -  Geometry  of  a  Coordinate  System  Fixed  to  a  Rotating  Blade 

In  this  geometry,  the  blade-reference  line  Is  any  reference  line  one  prefers. 

It  could  be  a  curve  through  the  midpoint  of  the  blade  chord  line  (curve)  on  the 
cylinder;  It  could  be  conveniently  the  curve  of  aerodynamic  centers  for  llftlng- 
llne  calculations;  or  It  could  be  the  leading  edge.  However,  the  choice  of  refer¬ 
ence  lines  should  be  judiciously  made  to  avoid  computational  difficulties.  Including 
the  Introduction  of  some  singularities  (e.g.,  the  leading  edge  of  conventional  pro¬ 
peller  blades  would  have  a  slope  of  Infinity  referred  to  chordwlse  distance). 


In  Figure  3,  defines  the  skew  angle  and  1^  defines  the  total  rake*  is 

the  pitch  angle  (the  orientation  of  nose-to-tall  helix  at  radius  r),  is  defined 

2  1  ' 
along  the  helix,  and  ^  is  perpendicular  to  ^  and  lies  in  the  cylinder  defined  by 

2 

r  =  a  constant.  Of  course,  ^  also  represents  a  helix  with  a  slope  that  Is  the 

negative  reciprocal  of  the  slope  of  at  the  point. 

1  2 

^  ,  r,  ^  Is  a  locally  orthogonal  set  but  does  not  represent  a  general  orthogo 

2  1 
nal  coordinate  system.  The  direction  of  ^  Is  in  the  blnormal  direction  of  the  ^ 

helix,  and  the  principal  normal  of  Is  in  the  negative  r  direction.  Even  though 

this  set  Is  not  a  general  orthogonal  curvilinear  coordinate  system,  one  still  can 

establish  the  equations  of  motion,  continuity,  etc.  In  this  coordinate  system.  The 

metrics,  Chrlstoffel  symbols,  etc.  Involve  rather  lengthy  algebraic  manipulations, 

'.'hlch  are  summarized  here. 

Blacle-Orlented  Coordinates: 


y^=  ^^cos^p  = 

y*=-r$in«* 

y*=  r  cosff 


ay^  .  . 


|.y-  =  i^+a0;  =  c 


=  -  COS^p 

where  0=  0^^  +  0,  +  +  ^  sin0p_  ^  +  “ 


1  2 

a  =  ^  costpp  sin0p 

(3  ~  ^'sin  0p  -  ^*cos^p 

and  1'  =  di  /dr;  0[  -d^, /df;  Op  =d0p/dr 


-  cos  6  cos  *Pp 


?)!’=  -  siaS -  r  cos  tf  [9;  -  •  -  ‘i"®  - ' » 


^-=  -coid  sin 0p 


^  =  -  sin0cos  (ftp 
6(' 


( 10  continued) 


costf  -  r  sintf  jej-  ^  =  cos®- r  B  sin® 

dr 


^  =  —  sin®  sin^p 

d{> 


3y"&y* 

From  equations  (10),  the  metric  summed  on  k  can  be  formed: 

ax' ax’ 


1  Csin*^p+  rBcos4>p  0 

Oji  =  Csin<^p+ rBco8<^p  l+r^6*+C*  “Cros<^p4  rB sin^p 


0 


—  Ccos^p  +  rBsin^p 


1 


where  C  =  i;+a^' 
and  h  --fj  - 


The  value  of  g  =Jg::|  —  1 


For  convenience,  some  combinations  occur  so  often  they  are  given  another  designation: 

Csin^p-t-rBcos^p  >K 

Ccos0p-  rBsin^p>  L 

When  these  are  Introduced,  the  metric  becomes: 


j=  K  HK*+L*  -L 

0  -L  1 


(lla) 


The  associated  metric  tensor  is  ■  G  '^/g* 

where  is  the  confactor  of  Because  g  “  1,  the  g^^  are  simply  the  cofactors 


They  are: 


g"=  l  +  C*sinVp+  r^B*cos*0p  +  r  BC$in20p 
g'i=g?t=  -  C  sin0p- r  Ecos0p 
gi3.g3i_j^B^^*jjn20p-  rBCco$20p 

g”  =  g”  =  Ccos^-  rBsin^p 

g'*  =  I 

g-’  =  J  f  Ccos*0p  -r  BCsin  20p 


Introducing  K  and  L  as  previously  defined,  one  obtains  the  tensor 


1+r 

-K 

-KL 

-K 

I 

L 

-KL 

L 

14  l" 

(12a 


Next,  one  must  calculate  the  Chrlstoffel  symbols  of  the  first  kind: 


(13) 


from  which  the  Chrlstoffel  symbols  that  appear  In  the  equations  are  derived • 


'il  =9'“[ii.l.] 


(14) 


(Note  that  T  Is  not  a  tensor;  some  authors  prefer  the  symbol  /  so  as  to  avoid 


any  confusion  from  the  so-called  Princeton  notation) • 


i:} 


[n,j]  =  o 

[11,3]  =[13,1]  =  [31,1]-  0 

[3  3, 1]=  [3  3,3]=  [13,3] -[31,  3]  »  0 

[i2.1]=[2M]  -  0 

[2  3.3]  =  [3  2,  3]  =0 

0'«2]  =  |^;*  =  |^[csin^P+rBcos^p]  = 
[3  3,2]  =  |S3”=  ^,[-Ccos^P+  rBsin^p]=  " 


ri3.2]=r3K2l  =  1  =--1. 

;  .  J  L  O  2  ^  2  r 

f)2.3]=f2i.  3]=J.  ^n_  =-4L 

2  hi'  hi* 

[32,»]  =  [23,i]=l  =^p 

[2  2  ,  =  _  jB+r^+^sin^p  -L<Ap-t-  Bcos^p 

3r  23^'  dr  dr  dr 

=  i *  sin  ♦(2«,'+ r«*  Jcos<^p+  ^.Vp  -  ^ '(♦p  )* 

[l3,  3]=[21,2]  =  1  |9;»i  =  -Bcos^p+^Ccos^p-rBsin^pj^p 

=  —  Bcos^p+  l^p 

=  cos^i-i^^pcos^p-  rj/^p  sln^p 
^cos*»p-i-4^in»pCos<»p  ^ 

^  ^*$in2»p  -  ^Vos2»p 

r 

[23,2]  =  [32,2]  =  -  B  sin0p  +  [C  sin0p  +  rBcos  0p]0p' 

=  ~B  $in0p  +  K^p 

=  -0,sin^p+l'^pSin^p4r5i0pCo$^P+£-£°*^Eiii! 


^p+(2^,'+  r^/j  sin^p-i 


K^p  f  B^in0p 


W  I 


^'^P  j^^^'cos«^p+  ^*in  <pp)-re';(^\\n<p 
^^coitppj  +-^j|^’ J  -  (f’^sin20p  -2^V*  os2^pj  j 
+  ^P  -r0;(^^in0p"^*cos<^p)  -^f(^^'‘--Os<^P+^*sin«^p)  + 

-J-  C1jLl(^_1  sin20p- fV*cos2<^pj| 

+(<^p)  -1t(^^*'''‘^p"^*<=os*^p)  ~  ^^i(^^'cos0p+f  sin<^»pj 

+  7[^(^’ )*"  (^*)^cos2<^P  +  2^V*sin20pjj 


The  Chrlstoffel  symbols  of  the  second  kind  are: 

r^1  ^  J^t  —  1^  $  I  n<j>pCO  S^D 
*  ^  3  I  f 

/rj=  /::  =  KBcos^p 


-Bcos0p 

1^1=  yr*  =  -  Bsin^p 

y^i_  pi  _  _  sin</vjCOS^ 

“  1-0P  +  Bcos0p  +  rKB* 

o  r 

/;'=-— -K0p  +  Bsin0p-  rLB* 

Sr 

r;^=-Li2^ 

pi__  COs’0o 
r 

77’= 

r 

p2_  sin’0p 
^33“  r 

/r,-/7'=  -0'p-LBcos0p 
/^j  /j  j  -  +  K  Bsin0p 


Looking  at  the  equations  of  continuity  and  motion,  (4)  and  (9a),  we  see  that  some 

derivatives  are  needed.  The  derivatives  of  the  Chrlstoffel  symbols  are  given  below 

12  13 

Derivatives  with  respect  to  ^  and  ^  or  x  and  x  coordinates: 


^  cos^0p 
S^«  S^'  r 

^ ^  $in<AKOs^ 


cos*0p 

_^]^sin0oCos0D 

be  ~w~be  r  — 
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?,Bcos^p+K|?fo,«p  ipl,  |^  =  |^Bcos0p+k||,co,« 
■ /bL  p  ■  1  ^  \  siprtp  3/7,  d/T,  dl  sind'oCOsd'p 


_  =:-^_sin^p£os^  ^^3  ^/7? _ /.^B+L^ 

“af'  a^'  f  ^ 

-  A/32  —  -  ( ^'•-B  +  L—  I  ^ _ ^  sin0p< 

a^'  \af’  afV  ^  5F  r 

.  sin*0p  3  =-^sin*‘^p 

a^'  df* 


^2  = 

a^' 


-_  ai  sin*0p 

a^' 

^^21  =  -  ^  COS<^p 

a  4'  a4 

=  —  —  sin<^p 

a^’  a^’ 

i/r!  =0 
a4^ 


^r*2  =  '^^1  =  -  ^5cOS<^p 

=  ^^  =  -|§sin0p 

d(*  af* 


1/7;  = 

^/72  = 

ar 

^^72  ... 

a^^ 

a/r? 

ai' 

a/;,^  _ 
a^' 

i£i. 
a(’  ■ 
a  773 
a^’ 

a/7f 
a«'  ' 

a/;; . 
a^‘ 


.  ^  0'+L®  cos0p  +  r  ^b’+sKESB 

'a^'ar  a^'  a^  a^'  a<' 

' 

:  a  K  a L  ^ ' -f- a B cos<iQ  +  ^  a  4  2k b^b 

b^dr  **  a^^  a^* 

P_ai^_  |K0'+^^sin0p-r(^B*+2BL^l 

a?  a^  *( 


-aT^- 


f 

|i,B>..Bt|B| 

tr-B’^^BLai) 


aK  sin^<^p 

a^’  ^ 


a/7’2  --/iko+L^^ 

^"a?^"  W  aW 

=^--0 

a^’ 


a/r7  _  aL 

ar  a? 

^=0 

a^’ 

ak"!  iH. 

i:iF" 

a/r’  a/rj  /aLg 
a^^  Va^» 

^  a«* 


•  «• 


<3/^1 


"V 


dK  \ 

i?/ 


£  2 

Derivatives  of  with  respect  to  r  (the  x  coordinate) 

|£l=  iKco^  -  K^^sinVos^p  ■«■  cos*<^  j 

_  d/;,'  aK  s\n4>pcos<hL  ^k(  cosVp-sin^<^p  a;:  sin«^pcos<^p> 
ar  ar  ar  r  \  r  r*  > 

a/?!  .  /^Q4.j^^\cos'#>p-KBsin'#>p'#'p 

Br  Sr  \ Sr  d '^  / 

_  /^p  4. 1  iJ?.\sin^p  -BLcos'/>p<#»p 

Sr  ar  Vsr  Sr  / 

^  sln^cos*^  4.  L  sin»/>p  cos«^p  4.1^/ sin*^p  ~  cos^^pWp 
df  "Sr  Sr  r  »■*  ' 

-  l^)sin^^p-  ^sirr^cos^p  i>p 
^  r  V  r^  r  /  ' 

2  7 

^Zlz=  =r  -  ^cos^pH-  Bsin'^'p  ^p 

9r  Sr  Sr 


^Ju  —  ^^3 ^sin*^p-Bcos<^p  ^p 

Sr  Sr  Sr 


2rB  — 
Sr  Sr 


Sr 

nil 

Sr 

ml 

S  r 

2 

nil—  2  cos«^psin*^BA'  ■  cos^*^p 
Sr  r  '*  P 


»SK 

Sr 


I  Dsin</>pVp-L<?^P  f fcos<#>p  +  2rKB)i  KB^-rS’- 
SrSr'VSr  /  Sr\  / 

=  (l  -J|l)cos’<^p^  iLcos^pSjn^p^p 
Vr^  r  Sr  /  r 
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/idK  KM 

|«in2<^p4-2K8(n^pCO_S^ 

V  rar 

bf\7  _ 

ML 

a  r 

ar 

ar?  = 

ay?? 

ar 

■v  1 

a  r 

^  I 

_ 

^^3  2 

=  +  j  sin<^p+KBcos'Pp'Pp 


n>;  following  are  given  to  express  L,  B,  K,  and  associated  derivatives  In  terms  of 
the  basic  geometry: 


B  =  SL  _ 

B  ^.77-7 


a  =  ('cos«^p+  ^*sin^p  fi  =  ^sirc^p-  ^*cos<^p 

L  =  (  +«‘^P )cos<^p- ( rd,'-  7-^^p)sin'/>p 

K  =  +a<t>p  )s\n<t>p  +  {rei-  7  -04>p)cos4>p 

^^4>'+  cos<^DSin*^p  ^  =  -  Sln^Pipp  -  cosJ>p 

f  c-V  r 

dL  _  sin^</>p  dB_  cos^^'_  sin«^n 

r  r  ^ 

dK _ cos'«/>p  ^  ^  g;  lA^P  +  ?«  -  9S!P'bS  ,;g<j 

d4'  r  dr  r 

^  <^3  -  cos<^pSin<^p 

dV  r 

—  -i.';sin<^p-fr<+  f  ".'\sin<^p-  ^sin<^p  1  ^'('/’p)  f  ^'<#>p  - K'/'p 


=  ^jsini>p  +^rel'+0i+-pjcos^p+  ^»*cos'^p-^’W)+4Vp  +  L^p 

=  sin'^p+^r5;"+2e,"  - -^)cos'^p+2'^p  i"cos'/>p - (rd^'+ei)sin^p 

_  !  ^  ']*(♦')[-  iTSln‘/»p  +  r0'cos^p+  Pos^p-^sin^p)j 

+  CVp"+  i^cos'^p-^rd^-  ~jsin<#»p 

0,=  t"’cos«#>p-^rC+20^^)sin.^>p+2<#»;  -i;sin«/>p- (re:+d.')cos^p 

+  .^cos<^p|  |_  i^cos*#‘p  +  r0^sin«#>p  -  ■^(o«sin*#>p +/ffcos<^p)| 

+  ^’<#>p  +<^p  -^,j,sin<^p  +  ^r©,' -  “^cos'^pj 


_  ^jnfpcosst'^,^'  _  §1q!^p 
a^’ar "'  '■  ” 

_  d,"  +  c£i-¥iZSiQ^pd>'  _  sin<^pcos«^p 

a^’ar 

a^K  _  2‘'^in<^pCps<^j,  A'  cp.^^ 

TJ:: —  r  r^ 

a^  or 

aV  _  _  cos^4sss]n^4\xf^  ^  sjn^pcps^p 

a^'ar  "  ''  *■" 


Sane  inslpht  about  the  behavior  of  derivatives  of  geometric  elements  is  given 
in  the  following: 

^  and  g  ^  are  dimensionless; 

Christoffel  symbols  of  the  second  kind  are  like  -(1/r),  and  differentiation 

2 

gives  terms  like  (1/r  ) ; 

,  3 

higher  derivatives  are  like  (1/r  )  etc; 

and  only  second  derivatives  end  up  as  derivatives  of  r  alone. 

These  properties  will  he  reflected  in  the  boundary-layer  approximation  if  this 

coordinate  system  is  used.  For  example,  if  one  were  to  consider  a  Flat  Plate  blade 

2 

operating  at  zero  angle  of  attack  (something  like  Meyne's  work  [7]),  i  would  be  the 


normal  to  a  blade  which  physically  would  have  to  lie  on  a  single  hellcoldal  sheet* 
Note  that  the  leading  edge  does  not  have  to  be  at  0  «  0  and  no  restrictions  are 
placed  on  the  blade  shape. 

2 

In  this  case,  one  boundary  Is  at  ^  *  0»  and  with  the  usual  boundary  layer 

2 

assumptions,  derivatives  of  the  dependent  variable  with  respect  to  ^  are  considered 

3  2 

to  Increase  the  magnitude  of  the  term*  That  is,  the  q  velocity  in  the  ^  direction 

1  2 

Is  considered  much  smaller  than  q  or  q  .  Thus,  to  retain  all  elements  of  the 

2 

continuity  equation  (1),  differentiation  with  respect  to  ^  must  increase  the 
magnitude.  With  this  conclusion,  only  the  term  Involving  d^^^^^wlll  be  of 
prime  Importance  In  the  viscous  terms.  However,  use  of  the  ccntravarlant  form  In 
the  pressure  term  with  a  non-orthogonal  set,  as  we  have,  gives  rise  to  some  com¬ 
plexity.  That  Is,  — ('/p)  g*^  will  have  (for  our  case)  three  pressure  terms 

■ '  all  equations.  Thus,  for  a  boundary  layer  similar  to  a  flat  plate,  one  would  be 
advised  to  use  the  covariant  form.  This  will  not  cause  any  real  problem  because  one 
will  undoubtedly  substitute  the  physical  veloclty^components  in  an  appropriate 
manner.  (Note:  It  Is  not  necessary,  though.) 

Now:  q“=.  v'V^ — .  where  v°  is  the  physical  quantity. 

V9aa 

then  q,  = 

There  are  some  problems  associated  with  the  previous  coordinate  system  when 
one  seeks  to  use  them  for  the  general  boundary-layer  problem  over  a  propeller  blade. 
One  difficulty  Is  to  establish  the  boundary  and  to  represent  the  normal  for  a 
general  warped  blade.  Hence,  the  physical  arguments  that  are  associated  with  the 
direction  of  differentiation  along  the  normal  are  somewhat  suspect.  However,  except 
in  the  area  of  the  leading  edge  of  the  blade,  for  many  blades  the  previous  formula¬ 
tion  will  be  adequate  over  most  of  the  blade  (this  would  be  true  for  most  thin 
blades  lying  close  to  a  helicoid). 

ORTHOGONAL  SURFACE  COORDINATES  (CONSTANT  FRACTION  OF  CHORD) 

Another  system  that  removes  the  restriction  of  the  blade  thickness,  camber, 
arbitrary  warp,  and  rake  can  be  established.  This  has  some  Improvements  but  is 
not  sufficiently  general  for  all  aspects  relative  to  the  flow  vield.  One  serious 
problem  Is  that  the  components  of  the  flow  outside  the  boundary  layer  must  also  be 


couched  In  this  system,  and  this  may  cause  considerable  difficulty  at  times. 
However,  It  Is  possible  that  this  system  may  still  have  advantages  outweighing  the 
above  problem.  One  distinct  advantage  Is  the  anticipated  ability  to  easily  treat 
the  leading  edge  problem  directly.  Including  the  stagnation  point  region.  The 
general  transformation  is  summarized  In  the  following.  Figure  4  illustrates  the 
necessary  elements  associated  with  this  system. 


Figure  4  -  Constant  Fraction  of  Chord  Blad^  Geometry 

If  S  (or  S/D  in  a  dimensionless  system)  be  the  position  vector  of  a  point  on  a 
surface,  then  dS  lies  in  the  surface.  If  dS  =3  — represents  the  vector  descrip¬ 
tion  of  dS  (which  lies  In  the  surface),  the  vectors  lie  In  the  surface. 

Furthermore,  because  we  are  dealing  with  a  surface,  only  two  surface  coordinates 
are  necessary  to  describe  the  surface.  Thus,  let 

S=§{Xp,Xc?  only 

~  y*^l  where  _bj  are  the  usual  cartesian  unit  vectors. 


We  will  denote  these  two  surface  coordinates  by 


Because  represents  a  chord- 


wise  parameter,  curves  described  by  setting  equal  to  a  constant  represent 
constant  chord  curves  on  the  surface. 
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•y-7 


Thus,  one  has  /P)  (jy  -t-  ^1.^1  ^jy  (15) 

SXr  "  aXc  ® 

where  Max.  must  be  a  vector  lying  along  the  curve  of  constant x  .  Generally, 


ds  ^u'^ds 


direction  cosine  of  a  surface  curve  where  ^1=!-  Isa  direction 

ds 


d V®  ^ 

parameter  of  some  sort;  furthermore,  it  is  uniquely  determined  once  — ^  and  — nre 

ds  du' 

known. 


Now 


3y  *  2 

,  i  ~^r~a  where  are  defined  as  before.  Selecting  u  as  the 
ou  ^ 


constant  chord  direction,  one  finds 


a,  r.  ^  ^  '  X  J  i  +  g,{  X„,  X  J  i  +h,(x^  ,  X,1  k  (16a) 


SX, 


Now  we  could  make  £2  a  unit  vector  by  appropriate  scaling,  but  there  seems  to  he 
little  gained  by  such  a  move;  it  seems  better  to  let  £2  6e  as  defined,  and  then 


a,*  83  *92  2  =-f2*+92*‘*’^2* 


Hence , 


IS2l=V922 


(16b) 


(17a) 


Also,  we  may  define  other  vectors. 


Thus , 


Eoni 


(17b) 
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and 


-  a 


5/nir 


Now,  because  MzbI  also  lies  In  the  surface,  a  normal  can  be  defined  as 

SXc 


-  a|%l.  ai^.sj34i)  .a,- 

3Xe  a  X„  a  Xe 


Again,  £2  Is  not  a  unit  vector. 


from  which,  according  to  our  above  definition,  f^,  and  *^3  are  the  components  of 
In  the  Cartesian  system. 


933  =  a 


,.a,=  alVpJ.a^ 

3xr  a Xg 


Finally,  we  may  take,  in  a  right-handed  manner. 


12*33=3  I 


an  orthogonal  vector  tangent  to  the  surface  (call  It  the  u  direction) 


l,-a,=  9„ 


|a|  %1  K  /  al%l ,  a{%l 


Thus  we  know  all  of  Che  elements  of  the  metric,  and  the  coordinate  system  Is  a 

general  curvilinear  coordinate  system.  The  proof  Is  quite  easy;  one  has  to  take 

12  3 

only  other  Inner  products  of  80*  The  U  ,  U  ,  and  U  were  chosen  In  the  manner 

described  previously  to  comply  with  a  right-handed  system  and  to  fit  the  usual  U^, 

2 

U  In  the  surface  so  that  one  can  follow  the  development  details  usual  to  Gaussian 

coordinates.^  It  might  be  well  to  note  here  that  the  development  to  this  point 

applies  to  the  suction  or  back  side  of  a  propeller  rotating  clockwise  about  the 

y^  =  X  trailing  axis.  On  the  face  (pressure)  side,  n  will  be  the  negative  of  the 

above,  will  have  the  same  description,  and  3^^“  x  £_  where  £-  refers  to  the 

pressure  side  and  Is  equal  to  -a^^  on  the  suction  side  as  described  previously. 

Thus  £  »  -£  X  £-  -  same  as  before.  However,  the  order  of 

i  Jp  ->p 

right-handedness  Is 

£j^-v  5.3 £2  1  when  substituting  In  the  equations. 

This  assures  that  is  still  toward  the  trailing  edge. 

One  remaining  preliminary  exists,  to  form  the  Christoff el  symbols  one  needs: 


^  ^  ^  ,  ^Oii  5Xc  Syl 

Su“  Sxr  Sy*  Sxe  Sy*  c>u® 


In  this  equation,  there  will  be  no  difficulty  (other  than  algebra)  In  determining 


are  the  ^<*  previously  defined.  However,  one  will  need  the 

deriva’  i'.cs  — ^  and  7  f  In  order  to  complete  the  process.  The  required  steps 

Sy*  oyi 

are  contained  In  the  following  section  along  with  specific  values  of  the  other 
terms. 


FORMULATION  OF  THE  GEOMETRY  FOR  SURFACE  COORDINATES 

The  blade  surface  Is  defined  by  the  vector  S/D  as: 


D“  (d  5(^e-0  5)sin^p-|cos^p^i  -  I'sinfli  +|cosdJi 
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G»  *  c/D(Xe  -O.5)sln0p- E/Dcos^p 


Clilll 

+^^(x.-o.5)sin<^p-  |j^^cos0p'^‘^p(St)(Xe-P5)cos^p  +E/bsin0pj 

^X||  W  R 

Wtih  a  slight  modification  of  Of  and  0  from  the  previous  work,  one  can  write:  • 

a-  96(Xe-0,5)  cos^p4-  E/Dsin^p 
9^(Xc-a5)sin<^p  -  E/Dcos^p 

\  J 

where  the  following  have  been  substituted,  as  in  Reference  3: 

c/q(x^(xc-0.5) 

^*=E/o(xp,xJ 

Other  useful  formulas  are: 

M  =0;-2«  -  ^  (xc  -0.5)cosfpt  ^(5^p)sin»p)  I 


where 


With  this  shorthand,  one  has: 


where 


=Gi- 

a,  Q„J,  +  F„eaf  FF„cos<^if  -Famine)  k 

(3  Xfl  2 

(24) 


eg=-COS0  sinOK  (see  references  1  and  3) 

e--sinO|  +  cos0k 


^(costfj+  sindK) 
SXe 


1  ^ 
2  3xc 


-F. 


"ax. 


dfi 


2  dxc 


(26) 


aA) 


e, 


Finally, 


Of 


course,  we  need  this  In  the  cartesian  Y  system; 
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+ 

+ 


33“ 


cos$ 

si  n0 


si  nd 

2  axe 


( 


"axc 


j 

I  ~ 

(28) 

(29) 


,»ow  we  can  express  the  metric  using  either  vector  form. 


(30) 


Lastly,  we  need 


anH  ^ 

ay'  ay' 


Now, 


x„- Iz/D)' 

x«  =  05f  I  cos'#'pv/ly/D)Vlz/br|^3r»  ’(-y/z)-0b  -0, 

x-i- 

-t-  -gJ^sln+p 
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and 


3tc/D)  _  djc/b}  ^  ^  (^t/D)^" 

3y!  dx„  3y'  '  By'  "  'ay'*  dy'  ~  By' 

thus 


3y'  “ 


BXn 

2y/D 

_2y/D 

aXq 

ay»" 

/y*+  z* 

‘■Dx/2 

at/a‘ 

3Xp 

72  ID 

2z/D 

axp 

ay* 

/y*4-  z* 

■  “D  x/2 

a(z/p 

1  sin<^p 

axe 

sin<^p 

ay 

■  D  c/D 

atx/Dl" 

(C/D) 

^=-2sind 


z2cosd 


(31) 


aiy^a 


^Xe  1 

a(y/a'“(c/L)) 


2</>p  sin0 (  a  s\n4*p-  ;?cos</>p)-^—  sin0cos'#>p 
-  cos^cos<#>p*-XnSind  cos</>p^,'  +  2(i^^D)sin  6  sin'/*p 

cos<^p|^sin<^p)sind 

d  Xp  prL) 


But 


asin<^p  ->3cos<^p  =  E/D 
acos«^p+^Sin«^p  =c/D(Xg-0.5) 

and  If  the  radial  variation  of  of f set  ware  simply  E  *  c(Xj^)Y(x  )*,  then 

d(E/D)  _  d(E/c-c/D)  _  E  d(c/P) 
dx^  dxp  c  dxR 

iliis _ = -l-r2sind(F_COS<^P+Gj5in<#'p)-COSecOS<^p'l  (31  continued) 

i(y/D)  c/dl  "  " 


*Thls  assumes  that  the  blade  section  geometry  is  not  changed  radially.  Many 

designs  are  such  that  E  ^  f.  .  ^  (x  )  +  -^(x  )  Y  (x  ) 

D  D  c^  R-'  c''  c^  -  c^  R''  c^ 

d  ’'c 

and  hence,  the  expanded  form  In  the  bracketed  terms  for  (Involving  a  and fi) 

must  be  used.  o(y/D} 


Finally , 


axc  _axc  _  1 

a(y*/D)  a(2/D)  c/o 


-2<#>pCOS^j^otcos^p->9sin<^pJ+^  cos®  cos^, 

$ 

-sinecos<^p-  2  [iT/D]cos©sin«/>p  — Xpd.^cos^oj 


_  1  “  2  cosdrFpCOS^^G^sin</»p1— sindcos<^p 

”(c/D)L  -*  J 


axc  1 

a(z/D)  (c/D) 


(31  continued) 


We  next  write  down  the  geometrical  quantities.  The  were  previously  written 
'•own  and  will  not  he  repeated.  The  following  is  a  direct  result  of  the  coordinate 
-system. 


Because 


9  =  9i  1  922933 

gii=.S!i=  J. 

9  g>i 

9^'=  -It  '  9^  —  — -4r 

y  -  -  y  ^  y  g53 


The  Christoffel  symbols  of  the  first  kind  are: 


2  ax‘ 


I" ''4  57 


All  others  are  zero. 

The  Christoffel  symbols  of  the  second  kind  are  given  by 


^i:  =  9'“ljk.a]  , 


Using  the  previous  results* 


r;=  s''M=9“^||: 

f7r9"N=-l^gf 

f7;-r:rg"tuii“  rg^'iif- 
Thc:  =^'|’3-''i49"|% 
[7k::  =g"iu2]=^g’^ 

n’|.g”;22.21=Lg«gl> 

r;fn:-g”I23.21=ig1^ 


f7rg”l11.2I=-lg’’P' 

f7:=g^11.3l=-^=^J' 

f7;=g“l22,3=-lg”|2^’ 

P3.21-2^’  a? 

f7’.=i;’=g“ii3.3i49’'^ 

(r^C.’=9’‘l23.3l=1^”^ 

[7;=g>=D3.3a=^”^ 


There  should  be  no  local  points  that  have  a  singularity  problem  for  a  finite 
thickness  airfoil  section,  and  the  arguments  used  previously  to  justify  only  con¬ 
sideration  of  the  highest  normal  derivative  In  the  viscous  term  should  still  he 
valid.  The  significant  single  problem  in  this  system  is  to  provide  the  boundary 
condition  at  the  outer  limit  of  the  boundary  layer.  The  surface  Itself  Is  the 
other  limit  with  X  going  from  0  to  1  and  X  from  the  bub  to  the  tip.  This  surface 

— ^  3 

will  be  given  by  S/D  and  represent  (for  our  purposes'  =0.  However,  the  outer 
boundary  condition  (at  »  or  the  6  botindary  layer  limit)  is  that  the  velocity  must 
approach  the  velocity  of  the  Inviscld  flow,  nr  at  least  approach  It  asymptotically. 
Furthermore,  the  quanti  t  i  r:-  2q  k  iP  and  ^?x(^xR)  must  be  expressed  In  this  coordinate 
system.  Now  these  vectors  are 


2qxl?  =2 


i 

q* 

-Q 


0  0 


=  2r-(q‘cos^pt  q’5in</'p)5,.q'f?cos<#>pe,-q'f?3n<^peJ 


*See  Appendix  B  for  details. 
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and  i2x(i2*R>-i2rf,  needs  to  be  put  Into  the  new  coordinate  system.  The  appropriate 
components  can  be  determined  by  taking  the  Inner  product  of  the  given  vectors  with 
a  unit  vector  In  the  appropriate  direction,  e.g.. 


q 


gg 

|gol 


physical  component  of  q  In  direction  of 


Using  the  1»  and  e  as  perhaps  the  simplest  description,  one  has 


as  physical  components  of  where  V8^2  ~  |— 2( 

The  ccntravarlant  (needed  first  In  any  case)  Is 

a  physical  component  of  q  In  a  direction 

’  ‘  Vg;. 

.•8., 

922 

Thus  every  vector  will  need  to  be  dotted  with  the  unit  vectors  ^  ^  ,  and  to 

9n  9^2  933 

convert  the  Invlscld  flow  field  vectors  to  the  components  needed  and  also  to  ex¬ 
press  the  Coriolis  and  centripetal  accelerations  In  the  equations  of  motion.  One 
should  be  careful  In  the  process  to  nondimenslonal Ize  the  quantities  using  P  as  the 
length,  because  all  geometry  was  done  this  way. 

Another  observation  may  be  made  that  may  be  helpful.  Presumably,  most  sections 
will  exhibit  some  stagnation  point  for  only  the  chordwlse  velocity  components  at 
some  point  other  than  =  0,  probably  on  the  pressure  side  of  the  hlade  for  the 
design  condition.  This  Is  usually  a  good  starting  point.  There  Is  a  good  chance 
that  with  the  favorable  pressure  gradient  that  the  flow  on  the  suction  side  may  be 
laminar  to  near  the  j olnt  of  minimum  pressure  hut  turbulent  elsewhere.  See  Ref¬ 
erence  8,  page  644.  It  may  be  possible  to  use  an  approximate  technique  similar  to 

9  5 

that  of  Wild  and  Nager. 
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CHORDWISE  IMBEDDED  ORTHOGONAL  SURFACE  COORDINATES  (See  Figure  5) 

In  this  coordinate  system,  one  coordinate  lies  on  a  constant  radius  cylinder 
and  is  tangent  to  the  blade  section.  As  before,  one  has 

5/D=  Gi-  ^sindi*^cos^K 

Now  sindk)  (24 

One  base  vector  is  tangent  to  the  blade  section.  Let  this  vector  he  C^; 


3Xe  3Xe“  3x/  3Xc 


anii  thus,  one  has 


2  9Xe 


2  oxc  3Xc  3xc'  J 


-t|Acos«  -(F,|&-(^os« 

-cos^l^  -sine 

OXc  3Xc 


3  a  fp  dff  0o(  "1  • 
3Xc  .  ”  3xc  ^"3xc 

/jL  fia.lJM.1  . 

\2  3xc  .3xc_ 


cos«)  j 


3x^  3Xc 


cose-^ 

9x, 


-h^  -rJ" 

x|  3Xc  OXc 
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and  In  the  cylindrical  polar  coordinate^. 


Also , 


If  = 


then 


a>t 


-ci— 

”3X(. 


f  -  +a_oi_ 

-3“  axc'  axc  -« 


-L  ^ 
2  ^ 


2 


(35b) 


1^ 

2ax^ 

.  M- 
■  3^: 


•cosd 


(35c) 


h  =4"'-sind 

3  aXc 

d)^  3Xb 

In  addition  to  the  above,  terms  like  »0yi  needed.  These  are  the  same  terms 

previously  defined.  The  form  of  the  Chrlstoffr 1  symbols  Is  again  the  same;  the 

difference  will  be  In  the  definition  of  and  g^^,  which  obviously  differ  from 

the  previously  defined  coordinate  system. 

Because  the  g.  ,  and  g^^  are  functions  of  x  and  x„ 

IJ  °  c  R 


ag.i  -dgjj  ax,  ,  agu  ax,  ay^ 
ax“  "axg  ay*  ^  ax*  ay  ax® 


The  dy^dx  ®  are  the  ay^au®  f^,g^  riven  on  page  34  for  this  coordinate 

system.  Thi  ax^  and  3x^/9  y*  are  given  In  the  previous  coordinate  system. 

The  ^  g.i .  g"  will  Involve  the  same  derivatives  in  rhe  present  system  (here) 

as  In  the  previous  surface  system.  Again,  because  a(S/0j/aX„and3(S/D)/3Xc  have 
the  same  directions  as  before,  n,  to  be  the  outward  normal,  will  be  the  negative  of 
equation  (29)  page  27  and  to  continue  to  have  c  the  same  (outward),  we  will  have 
to  reorder  t^ e  coordinates  for  the  pressure  side  as  follows; 

np=  -  ( back  or  suction  side] 


34 


80  that  Is  defined  as  before.  The  right-handed  system  will  now  require  that  the 
order  is 

fip— ^  C,  — ^  £,  pressure  side 

and  — *"  suction  side 

and  Hp  =  rip  =-\  =-ng 

The  appearance  of  the  coordinate  system  is  shown  in  Figure  5.  The  boundary 
condition  at  the  outer  edge  of  the  boundary  layer  here  may  be  very  close  to  the 
specif Ictlon  that  the  direction  is  a  streamline.  If  this  is  true,  there  may  be 
some  simplification. 


Vigure  5  -  Chordwlse  Imbedded  Coordinates 

The  equations  are  now  complete.  Appropriate  geometric  configurations  can  now  be 
selected  and  Individual  terms  examined  to  derive  a  set  of  equations  appropriate  for 
the  boundary  layer  flow. 

FINAL  REMARKS 

One  should  be  able  to  perform  houndary-layer  calculations  with  the  coordinate 
systems  developed.  A  zero  pressure  gradient  solution  may  he  made  for  a  skewed 


blade  similar  to  Meyne*  However,  the  real  blade  Is,  of  course,  the  necessary 

problem*  The  leading-edge  region  Involving  stagnation  point  flow  In  the  u^ 

direction  (perpendicular  to  the  constant  chord  line)  probably  should  be  approached 

first*  One  usually  can  expect  some  uncoupling  of  the  equations  In  this  region,  and 

a  laminar  solution  may  be  both  amenable  and  practical  because  there  are  favorable 

pressure  gradients  In  the  vicinity  of  the  stagnation  point,  which  may  be  particularly 

strong  on  the  suction  side*  Whether  a  series  or  a  polynomial  approximation  Is  used 

may  not  make  too  much  difference*  The  ensuing  laminar  solution  on  the  rest  of  the 

blade  may  be  academic  and  of  little  practical  use,  but  some  model  for  the  transition 

region  and  turbulent  region  beyond  can  be  applied* 

In  the  case  of  turbulence,  one  must  add  the  Reynolds  stress  terms  that  appear 

as  a  result  of  the  momentum  terms  In  equations  (9a)  and  (9b)*  Because  q^  and  q^  do 

not  appear  anywhere  else  as  nonlinear  terms,  this  Is  the  only  addition  necessary* 

The  choice  of  coordinates  for  the  turbulent  boundary  layer  problem  may  be  the  latter 

surface  coordinate  set  ( x  >•  a  constant)  because  physical  evidence  seems  to  Indicate 

R 

the  flow  Is  essentially  In  that  direction*  It  may  be  that  the  original  coordinate 
set  may  be  used  with  the  appropriate  pressure  gradient  included*  This  would  be  con¬ 
venient  because  no  component  construction  would  be  needed  for  either  the  pressure 
gradient  or  the  matching  fluid  dynamic  velocities  as  one  proceeds  to  the  limit  of 
the  boundary  layer*  The  appropriateness  of  this  latter  possibility  may  have  to  be 
determined  by  two  calculations  using  the  surface  coordinates  in  one  computation  and 
the  references  surface  coordinates  In  the  other* 
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so  chat  is  defined  as  before*  The  right-handed  system  will  now  require  that  Che 
order  Is 

— ^£3  pressure  side 

and  — ►•£3 — ^  £,  suction  side 

and  Op  =  rip  =-03  =-n3 

The  appearance  of  the  coordinate  system  Is  shown  In  Figure  5.  The  boundary 
condition  at  the  outer  edge  of  the  boundary  layer  here  may  be  very  close  to  the 
speclflctlon  that  the  direction  Is  a  streamline*  If  this  Is  true,  there  may  be 
' ome  simplification* 


Vigure  5  -  Chordwlse  Imbedded  Coordinates 

The  equations  are  now  complete.  Appropriate  geometric  configurations  can  now  be 
selected  and  Individual  terms  examined  to  derive  a  set  of  equations  appropriate  for 
the  boundary  layer  flow. 

FINAL  REMARKS 

One  should  be  able  to  perform  boundary-layer  calculations  with  the  coordinate 
systems  developed.  A  zero  pressure  gradient  solution  may  be  made  for  a  skewed 
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blade  similar  to  Meyne.  However,  the  real  blade  Is,  of  course,  the  necessary 
problem*  The  leadlng~edge  region  Involving  stagnation  point  flow  In  the 
direction  (perpendicular  to  the  constant  chord  line)  probably  should  be  approached 
first.  One  usually  can  expect  some  uncoupling  of  the  equations  In  this  region,  and 
a  laminar  solution  may  be  both  amenable  and  practical  because  there  are  favorable 
pressure  gradients  In  the  vicinity  of  the  stagnation  point,  which  may  be  particularly 
strong  on  the  suction  side.  Whether  a  series  or  a  polynomial  approximation  Is  used 
may  not  make  too  much  difference.  The  ensuing  laminar  solution  on  the  rest  of  the 
blade  may  be  academic  and  of  little  practical  use,  but  some  model  for  the  transition 
region  and  turbulent  region  beyond  can  be  applied. 


In  the  case  of  turbulence,  one  must  add  the  Reynolds  stress  terms  that  appear 
as  a  result  of  the  momentum  terms  in  equations  (9a)  and  (9b).  Because  q^  and  q^  do 
not  appear  anywhere  else  as  nonlinear  terms,  this  Is  the  only  addition  necessary. 


The  choice  of  coordinates  for  the  turbulent  boundary  layer  problem  may  be  the  latter 


Surface  coordinate  set  (x  =  a  constant)  because  physical  evidence  seems  to  Indicate 

R 

the  flow  is  essentially  In  that  direction.  It  may  be  that  the  original  coordinate 
set  may  be  used  with  the  appropriate  pressure  gradient  Included.  This  would  be  con¬ 


venient  because  no  component  construction  would  be  needed  for  either  the  pressure 


gradient  or  the  matching  fluid  dynamic  velocities  as  one  proceeds  to  the  limit  of 


the  boundary  layer.  The  appropriateness  of  this  latter  possibility  may  have  to  be 
determined  by  two  calculations  using  the  surface  coordinates  in  one  computation  and 
the  references  surface  coordinates  In  the  other. 
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Appendix  A 

Some  Tensor  incepts  and  Identities 


We  define  a  tensor  by  the  way  the  quantity  transforms,  for  example,  B(i,x)=3y*/dx‘  A(o,y)  covariant 
transformation.  Sometimes  B(i,x)»dx'/dy*A(ar,y)  contravariant  transformation. 

This  does  not  make  much  sense  unless  one  looks  at  some  examples.  Obviously,  a  scalar,  f  say,  is  the  same 
in  any  system,  but  look  at  the  derivative.  3f/3x‘ =3f/3y*‘9yV3x'  This  is  a  covariant  transformation.  We 
use  F,  to  represent  a  first  rank  tensor  that  transforms  in  a  covariant  manner  and  F'  for  a  first  rank  tensor 
that  transforms  in  a  contravariant  manner.  As  one  might  expect,  the  df/dx'  is  in  some  way  related  to  a 
vector  in  the  x'  direction,  which  is  true.  The  precise  meaning  of  F^ ,  F'  in  terms  of  the  usual  vector  com¬ 
ponents  will  be  discussed  shortly. 

First,  however,  one  must  understand  that  higher  rank  tensors  than  the  first  exist.  Typically,  in  fluids  one 
may  get  to  third  rank  tensors,  in  solids  all  the  way  to  fourth.  Tensors  may  be  mixed,  e.g.,  R]^|(y)  is  a  fourth 
ranked  tensor  and  transforms 


ay' 

ax* 


ay'  ay**  ay*  ' 


N'ote  that  the  transformation  is  by  products  of  first  partial  derivatives.  Furthermore,  one  should  sum  on 
the  repeated  sub  and  superscript  and  a  superscript  in  the  denominator  can  be  considered  as  a  subscript  in  the 
numerator.  Coordinates  are  labeled  with  superscripts  because  they  transform  in  a  contravariant  manner. 


Differentiation 


If  we  differentiate  a  f/  a  x '  above,  we  get  , 

aM  ^  an  ^  ^  ^ 
ax'ax'  ay*ay'  ax'  ax*  ay*  ax'ax* 

and  hence  the  expression  does  not  transform  as  a  tensor.  To  overcome  this  problem,  proceed  as  follows: 

aB,(y)  ax*  ax'  aA  a*x*  . 

— lill  =  -  -  — 1  +  -  A 

0 

ay'  ay'  ay*  ax'  ay'ay' 

If  we  look  at  the  Christoffel  symbols,  one  can  devise  the  following: 

a^  X*  _  d)r  ax' 

ay'ay*  ”  '*  ay**  *  ^  ay'  ay' 

Inserting  this  in  the  top  equation  gives 


3B,(y) 


-  r‘B 

y  *i  •« 


dX’  ax»  /  3^  T'  A  ^ 

a7‘  fly'  " «  •#  >) 


which  transforms  as  a  tensor. 


We  designate  .^  =  3 / 3y'  -  PJ  B^  as  covariant  differentiation. 

l£the  vector  was  a  contravariant  component,  then  B',,=  3B'/3y*+  PJi^  B^  Let  the  physical  component 
of  A  be  Ap  and  be  expressed  in  contravariant  form  by  Aj,  =  \/~^  A'  (not  summed  on  i).  The  contravariant 
components  of  A  are  A'  =  Aj/vAg^  (not  summed  on  ‘ )  where _A‘ ,  the  physical  components,  represent  the 
sides  of  a  parallelopiped,  of  which  A  is  the  diagonal.  Thus  if  q  (U,  ,  V^),  where  U,  ,  and  V^are  the 

physical  components,  then 
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in  cylindrical  coordinates. 


The  covariant  component  of  A'is  obtained  in  a  similar  manner: 

A,-g„  A'  =Vg;rA'p  (not  summed  on  i) 

The  curl  of  a  vector  Q  represented  by  V  x  Q  has  components  as  follows.  (Note:  we  use  e-system  here.)* 


♦e  '"  is  1  i  #  j  #  k  #  i  and  i  -*  j  -*  k  is  cyclic  in  order,  e.g.,  123 
=  — 1  i  /  j  #  k  #  i  and  i  -*  j  -►  k  is  acyclic  in  order 
=  0  if  i  =  j  or  j  =  k  or  i  =  k 
e  is  the  same. 


s 


where 


q 


I  ’k 


-  r*  q 

)k 


The  vector  or  cross  product  of  two  vectors  is 

H  =  Q  X  q 


and  has  components  H'  =  0'"'  12 
or  H  =  6  Q**"* 

Of  course,  H '  =  Q ''  H ^ ,  so  if  contravariant  forms  of  the  vector  are  known,  one  may  write  H  ^  instead  of  H 
Among  other  useful  things  concerning  the  use  of  the  0  system  is: 


One  literally  has  the  choice  of  using  covariant  or  contravariant  components  to  determine  the  products 
formed  by  various  vectors.  The  only  reason  for  using  contravariant  forms  is  if  one  works  with  the  equations 
of  motion  and  continuity  in  this  form;  then  it  is  very  easy  to  convert  to  the  physical  form.  Neither  form  is 
really  more  appropriate  than  the  other. 

Often  the  equations  are  manipulated  in  one  form  and  then  converted  to  physical  components  before 
actually  solving  the  equations.  This  is  perhaps  most  common,  but  not  necessarily  most  efficient. 

Continuing,  we  have  a  triple  vector  product  involved  in  our  equations  because  of  the  rotating  coordinate 
system.  (Note  that  we  have  introduced  the  general  Kroenecker  delta  )! 

Q  X  (Q  X  R)  for  which 

K'  =  d"Q 

pq  I 

In  the  X,  r,  6  system, 


-  Q  =  Q’  only,  because  9,,  =  1 


and  -  Q  =  S2j 

R  is  the  general  position  vector  and,  because  of  the  rotation  only  about  X,  only  the  ©^component  of  R  is 
important,  although  both  R^  and  R^  (but  not  R^ )  exist. 


=0  unless 


i  /  j  sfc  I  i 
p  #  q  #  r  #  p 


=  ±  1  i,  j,  k  are  repeated  in  any  order  in  p,  q,  and  r.  If  the  number  of  permutations  (one  at  a  time)  is 
even  to  put  pqr  in  the  same  order  as  ijl,  the  plus  is  used;  if  this  number  of  permutations  is  odd,  the  minus  is 
used. 


(Note:  is  also  1) 


R  =  R>‘e.+  R"e, 

K'  =  d''  fi  Q' R"  -  R' d'’=  0 
1||  ^  11 

K*  =  Q,  Q’  R”  q  =  2  and  one  permutation  is  needed  to 

”  get  symbols  in  the  same  order. 

R*  =  -  r 

For  further  information  one  should  refer  to  [6]  or  [10]. 
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Appendix  B 

Christoffel  Symbols  for  Blade  Surface  Coordinates 


Some  of  the  details  to  carry  out  the  determination  of  the  Christoffel  symbols  in  the  blade  surface  coor¬ 
dinate  systems  are  carried  out  here.  As  shown  in  equation  (22)  repeated  here, 

££ii  =  i£ii  (22) 

du"  3*^  ay'  au*  a*_.  ay'  au* 

It  is  necessary  to  differentiate  Fp.Gp,  ao/ax^,  and  a/J/ax^  which  are  contained  in  the  g,,  (equations 
(30)).  The  aXp/Qy^and  3 x^ /ay  are  given  in  equations  (31)  and  the  ay' /au*are  the  i,  j,  and  k  com¬ 
ponents  in  equations  (24),  (28),  and  (29);  that  is. 


ay’  .  ay*  .  ay* 

a  =  b  ,  =  I  =  I  + 


—  •  —  J 


au*  -  au*  ^  au* 

=  (oi  Qoi  (’ak 


(Bl) 


Equatio.is  (B2)  and  (B3),  combined  with  (30)  and  (31),  will  allow  one  to  complete  the  determination  of  the 
Christoffei  symbols,  and  thus  the  algebraic  form  of  the  equation  to  be  solved. 
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Appendix  C 

Cylindrical  Polar  Coordinates 


A  common  appearing  coordinate  system  is  given  in  the  literature  and  listed  below  for  completeness. 


1  0 
0  1 
0  0 


;  g  = 


9ii 


“1  0  0* 
0  1  0 

.0  0  ^ 


The  only  nonzero  Christoffel  symbols  of  the  first 
kind  are: 


[32,2]  =  [23,2]  =  -[33,2] 

=  r 

The  Christoffel  symbols  of  the  second  kind  are: 

r;  =  g  ^”[ij,a] 

n,=  n3=g“[23,3]=[1/r^]r  =  1/r 
r^=g“[33,2]  =  -r 


All  others  are  zero. 
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OTNSROC  ISSUES  THREE  TYPES  Of  REPORTS 

1.  OTNSROC  REPORTS,  A  FORMAL  SERIES.  CONTAIN  INFORMATION  OF  PERMANENT  TECH 
NICAL  VALUE.  THEY  CARRY  A  CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 

2.  DEPARTMENTAL  REPORTS.  A  SEMIFORMAL  SERIES.  CONTAIN  INFORMATION  OF  A  PRELIM 
INARY,  TEMPORARY.  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE 
THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3.  TECHNICAL  MEMORANDA,  AN  INFORMAL  SERIES,  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INIEREST.  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN 
TERNAL  USE.  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  CODE  OF  THF  ORIGINATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CASE  BY  CASE 
BASIS. 


END 
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